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Abstract 

Friction Stir Spot Welding involves complex physical phenomena, which are very difficult to probe experimentally. In this 
regard, the numerical simulation may play a key role to gain insight into this complex thermo-mechanical process. It is 
often used to mimic specific experimental conditions to forecast outputs that may be substantial to analyse and elucidate 
the mechanisms behind the Friction Stir Spot Welding process. This welding technique uses frictional heat generated by 
a rotating tool to join materials. The heat generation mechanisms are governed by a combination of sliding and sticking con- 
tact conditions. In the numerical simulation, these contact conditions are thoroughly dependent on the used friction model. 
Hence, a successful prediction of the process relies on the appropriate selection of the contact model and parameters. This 
work aims to identify the pros and cons of different friction models in modelling combined sliding-sticking conditions. A 
three-dimensional coupled thermo-mechanical FE model, based on a Coupled Eulerian-Lagrangian formulation, was devel- 
oped. Different friction models are adopted to simulate the Friction Stir Spot Welding of the AA6082-T6 aluminium alloy. 
For these friction models, the temperature evolution, the heat generation, and the plastic deformation were analysed and 
compared with experimental results. It was realized that numerical analysis of Friction Stir Spot Welding can be effective and 
reliable as long as the interfacial friction characteristics are properly modelled. This approach may be used to guide the con- 
tact modelling strategy for the simulation of the Friction Stir Spot Welding process and its derivatives. 
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physics behind the process is still not thoroughly under- 
stood, since it involves coupled thermo-mechanical 
phenomena. 


Introduction 


Friction Stir Spot Welding (FSSW) is a solid-state spot 
welding technique derived from the linear Friction Stir 
Welding (FSW) process. This joining technique was 
developed and implemented by Mazda Motor 


‘Laboratoire de Génie Mécanique, Ecole Nationale d’Ingénieurs de 


Corporation in 1993 as an alternative to single-point 
joining methods such as Resistance Spot Welding 
(RSW).' Currently, despite being the most widely used 
welding process for vehicle structural assemblies, RSW 
may imply metallurgical issues associated with solidifica- 
tion for lightweight structural materials, like aluminium 
alloys.” FSSW process has received particular attention 
from both scientific and industrial communities. In this 
welding technique, the welds are created by plunging a 
rotating tool, into two overlapping plates, which are sof- 
tened by frictional heating and joined by plastic deform- 
ation.> Due to plastic deformation, FSSW has the 
potential to produce high-strength joints.* However, the 
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Scarce published research works were addressing the 
analysis of the heat generation mechanisms and their rela- 
tionship with the temperature evolution during welding. 
The analysis of heat generation was based on the conven- 
tional Coulomb law to describe the friction phenomenon 
at the tool-workpiece interface. The main reason for the 
scarcity of reports on such a subject is the complexity of 
the involved phenomena, which are very difficult to be 
experimentally probed. However, owing to the improve- 
ment of the computer power, researchers have developed 
numerical models to analyse the heat generation in FSSW. 
Finite Element Method (FEM) has been used for the 
numerical simulation of FSSW since about two decades 
ago. Several numerical analyses have been carried out in 
this respect and a wide range of numerical models were 
developed. 

Muci-Kiichler et al.° developed a 3D numerical model 
to simulate the isothermal plunge phase of the FSSW 
process. They focussed on the analysis of the flow stresses 
and the plastic deformation during welding. However, no 
thermal phenomena were analysed. Afterward, miscella- 
neous coupled thermo-mechanical models were devel- 
oped to gain more insight into the fundamentals of the 
FSSW process.° * Those models were used to predict 
the stress fields, the plastic deformation, and the tempera- 
ture distribution during welding. They were based on 
some simplifications that limited an efficient analysis of 
such a complex FSSW process. The comparison 
between the numerical models outputs and experimental 
results, such as temperature, torque, force, etc., was pre- 
sented without any analysis of the thermo-mechanical 
phenomena involved during the welding process. 
Therefore, the reliability of numerical modelling is still 
a debateable topic. 

Based on the abovementioned literature survey, it is 
likely to state that a successful FSSW numerical simula- 
tion could be established based on an appropriate selec- 
tion of the FE formulation, of a reliable constitutive 
material model, and of a friction model. In literature, 
special attention has been paid to the selection of the 
appropriate numerical formulation.” '* The combination 
of Lagrangian and Eulerian formulations was demon- 
strated to be the most effective for modelling the friction 
stir welding (FSW) process and its variants. On the other 
hand, the selection of the most appropriate material model 
for FSW is still limited. Generally, standard empirical 
models, such as Johnson Cook*!*:!4 and Norton-Hoff!~ 
'7 constitutive models, are frequently used. Besides, the 
conventional Coulomb’s law was used in the majority 
of FSSW numerical simulations. This friction law 
includes one coefficient of friction (COF) which is 
assumed as a constant value along the entire contact inter- 
face. Several values of COF were proposed for 
steel-aluminium pairs. Frigaard et al.’* and D’Urso 
et al.'? considered a COF of 0.4 as an average value 
from results in literature that reported on the FSW of 
AA6082-T6 aluminium alloy. Yang et al.”? used a COF 
of 0.577 in the numerical simulation of the FSSW of the 
AA2198-T8 aluminium alloy; but, no argument was 


given on how this value was determined. Likewise, 
Bagheri et al.'* assumed a COF of 0.35 for modelling 
the friction during the FSSW of AA5083 aluminium 
alloy. 

In the abovementioned works, the evolution of tem- 
perature and the material flow during FSSW process 
were presented. However, none of them addressed the 
influence of the contact at the tool-workpiece interface 
on the output results. The contact conditions at the tool- 
workpiece interface are dependent on the local tempera- 
ture and stress fields. Recently, Sanjeev et al.'’ and 
Chen et al.8 established a modified version of 
Coulomb’s law in which the limit shear stress of the 
welded material was introduced to control the stick-slip 
behaviour. In these works, the axial force, as well as the 
thermal history near the contact zone, were assessed. 
The obtained numerical results showed a good agreement 
with experimental ones. Scarce studies have used 
temperature-dependent COF for modelling the contact 
during welding. Khosa et al.*! employed a temperature- 
dependent friction coefficient determined using an 
inverse approach. They developed a 2D numerical 
model to simulate the FSSW of AA6082-T6 aluminium 
alloy and a comparison with experimental results was per- 
formed. They reported that the thermal prediction was 
satisfactory. On the other hand, Meyghani et al.”* devel- 
oped a mathematical formulation to calculate the 
temperature-dependent friction coefficient during the 
FSW of AA6061-T6 aluminium alloy. This formulation 
was assessed with various rotational speeds and showed 
its capability to predict the peak temperature. 

To date, only few analyses were achieved on the per- 
formance of different friction models on the simulation 
results of the FSSW process. Zhang et al.?* compared 
the classical Coulomb’s law and its modified version in 
modelling FSW of AA6061-T6 aluminium alloy. The 
authors considered a limit shear stress at the contact inter- 
face between the tool and the workpiece. They found for 
high rotational speed that the modified Coulomb law was 
more efficient in the FSW simulation. Dialami et al.** 
developed a modified version of Norton’s law implemen- 
ted in COMET”? a dedicated FSW software. They took 
into account the sliding velocity and the pressure distribu- 
tion at the tool-base material interface for modelling the 
contact conditions in the FSW of AA6063-T6 aluminium 
alloy. They compared Norton’s law and its modified 
version predictions and they showed that the modified 
Norton’s law has a great potential to reproduce experi- 
mental results. However, it is worth mentioning that in 
such cases, usually the agreement between numerical 
and experimental results was based on trial and error para- 
meter fitting of the model. Therefore, the study of the 
effect of friction models on the outputs predicted by 
numerical simulation was not dealt with thoroughly and 
it is of great necessity. The lack of reports in this field con- 
sists the main motivation of the current work. 

This work puts forward a fundamental basis for under- 
standing the effect of friction models on temperature evo- 
lution and heat generation mechanisms during the FSSW 
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of the AA6082-T6 aluminium alloy. A three-dimensional 
(3D) coupled thermo-mechanical model using a coupled 
Eulerian-Lagrangian (CEL) approach was developed 
under ABAQUS/Explicit software. For a large range of 
temperatures, experimental stress-strain curves were 
used to describe the base material behaviour during 
FSSW process. Three friction models were adopted, 
namely: the conventional Coulomb friction model 
(CLB), the modified Coulomb model (MDF-CLB), and 
the temperature-dependent friction coefficient-based 
model (COF(T)) in the numerical simulation of FSSW. 
These friction models were chosen to study their effect 
on the heat generation mechanisms in the FE analysis of 
FSSW process. Particularly, their effect on the tempera- 
ture evolution and both frictional and plastic energy dissi- 
pation was discussed. The findings of this investigation 
are helpful for directing the selection of the most appropri- 
ate friction model used for the FSSW simulation. The vali- 
dation of the developed model was carried out by 
comparing predicted results with those obtained 
experimentally. 


Thermo-mechanical modelling 


Geometry, mesh, and numerical formulation 


The numerical simulations were performed using 
ABAQUS/Explicit Finite Element code, which is well- 
established for handling non-linear problems. A fully 
coupled thermo-mechanical FE model was developed 
using the Coupled Eulerian-Lagrangian (CEL) approach. 


CEL formulation enables to selectively mesh the different 
model components, i.e. parts of the bodies undergoing 
large deformations are modelled using the Eulerian for- 
mulation, while the remaining are modelled using the con- 
ventional Lagrangian formulation. 

Figure | shows the built CEL model, i.e. the partitions 
(a), the dimensions (b), and the mesh topology (c). For 
saving computation time, only the central part (100 x 
100 x 10 mm) of the experimental samples was modelled. 
Moreover, to accurately capture differences in heat 
affected zones, a single rectangular sheet with a thickness 
of 10 mm substituted the two plates lap-welded. A similar 
procedure was used in the experimental tests. The idea of 
substituting the two lap plates with a single plate was 
based on Costa et al.”° research work. It was stated that 
since the two plates are in very tight contact during lap 
joining, the heat distribution is uniform through the two 
plates thickness. A 10 mm diameter pinless cylindrical 
tool was used in this investigation. 

The pinless tool was modelled as a rigid Lagrangian 
body and meshed with 13,271 four-node thermally 
coupled elements (C3D4T). The workpiece was modelled 
as an Eulerian domain, meshed with 58,240 multi- 
material thermally coupled 8-node Eulerian bricks. A 
refined mesh was assigned to the workpiece through the 
thickness direction (4mm) to enhance the accuracy of 
the computation near the contact interface. As shown in 
Figure 1, a fictious layer was placed on the top of the 
Eulerian domain, which was created for enabling the 
visualization of any material motion outside the work- 
piece upper face. No constitutive properties were assigned 


(a) | 


— 


Lagrangian 
Domain 
(C3D4T) 


Eulerian 
Domain 
(EC3D8RT) 


Figure I. The model description: (a) CEL model partitions; (b) part dimensions; (c) mesh topology. 
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to this top layer. The size of this domain was chosen to 
prevent any outside material flow, which would become 
invisible during the ABAQUS post-processing. The mate- 
rial flow visualization was based on the output variable 
Eulerian Volume Fraction (EVF).*’ An EVF equal to 
one indicates that the mesh element is completely filled 
with the workpiece material, while an EVF of zero indi- 
cates that the element is empty (void). The computational 
time required for performing each analysis is about 69 h 
on a 2.80 GHz Intel® Xeon® processor. 


Material behaviour 


In order to accurately describe the material behaviour in 
FSSW, the material model should take into account the 
temperature and strain-rate dependency. In the present 
research work, the AA6082-T6 aluminium alloy, a heat- 
treatable (precipitation-hardened) aluminium alloy, was 
used as base material. This alloy is frequently used in 
welding construction and well-characterised in previous 
research works from current co-authors. Typically, the 
mechanical properties of the AA6082-T6 aluminium 
alloy exhibit a significant decrease during welding, due 
to the dissolution of the strengthening precipitates (B 
family).’* In fact, starting at temperature 100 °C, the pre- 
cipitates exhibit further growth (overaging) that weakens 
the precipitate-dislocation interaction, which induces the 
decrease of the yield strength. This material behaviour 
should be considered in the FSSW process simulation. 
Hence, a thermo-elastic-plastic material behaviour was 
adopted for the current simulations. Stress-strain data 
points, for different temperature values, were introduced 
in the property module of ABAQUS/Explicit. The input 
data were extracted from the stress-strain curves presented 
in Figure 2. These curves were obtained from the experi- 
mental results reported by Leit&o et al.*? and Summers 
et al.2° According to Chen et al.* the strain rate has a 
minor impact on the flow stress and fracture behaviour 
of this AA6082-T6 aluminium alloy; therefore, its effect 
was neglected in the material modelling. 


Contact modelling 


The main objective of this paper is to use three friction 
models in the FSSW simulation and to assess their 


True stress (MPa) 


0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 


True strain 


Figure 2. Base material behaviour at different temperatures. 


ability in the description of the inherent thermo- 
mechanisms during welding. These models are as follows: 


e The conventional Coulomb model (CLB) which 
assumes that the critical shear stress at the tool- 
workpiece interface is proportional to the applied 
normal stress,*’ and expressed as 


Tcrit = PON (1) 


where 7,,; 18 the frictional shear stress, oy is the normal 
stress and yp is the COF. This coefficient is constant 
along the contact surface during the simulation.** In 
this work, a constant COF value of 0.3 was selected 
according to different works on the modelling of 
FSW of 6xxx aluminium alloys using the classical 
Coulomb law.??3° 

The modified Coulomb model (MDF-CLB), in which 
two different critical shear stresses may be defined: 
the first one is given by equation (1), which is active 
when the contact shear stress is still lower than the 
yield shear stress of the material (tyjeiq), the second 
one is given by 


Terit = Tyield (2) 


which is active when the contact shear stress reaches 
the yield shear stress of the base material. The yield 
shear stress of the material is estimated to be 


Tyied = (Gyiea | V3), (3) 


where Oyie is the base material initial yield stress. 
Using ABAQUS© software, only a constant value of 
the critical shear stress could be introduced in the 
Interaction module. Therefore, in the current analysis, 
the tyie adopted for the MDF-CLB model was an 
average value calculated using all the values of 
(Cyieia) presented in Figure 2. Hence, the critical shear 
stress is expressed as follows: 


Terit = MiN( HON, Tyieid)- (4) 


The temperature-dependent friction coefficient-based 
model (COF(T)). According to this model, the critical 
shear stress is defined as 


Terit = MT)on, (5) 


where p(T) is the temperature-dependent COF. The values 
of (T) used in the current analysis are given in Table 1.7 


Boundary conditions 


In the current model, both thermal and mechanical bound- 
ary conditions are introduced. The thermal boundary con- 
ditions were assigned by prescribing an initial ambient 
temperature to both the tool and the workpiece. The 
heat transfer between the workpiece and the surrounding 
air was simplified using a convection coefficient of 
30 W/m?C. Thermal conductance was set to the contact 
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Table |. Evolution of friction coefficient with temperature.” 


Coefficient of friction Temperature (°C) 


0.25 25 

0.23 50 

0.20 100 
0.19 150 
0.18 200 
0.15 250 
0.10 300 
0.05 350 
0.025 400 
0.015 450 
0.010 500 


Table 2. Thermal and mechanical properties of the workpiece 
and the tool materials. 


AA6082-T6 H13 steel 
(workpiece) (tool) 
Young’s modulus; E 70 210 
(GPa) 
Poisson’s ratio; v 0.33 0.1 
Yield strength (MPa) 325 1650 
Ultimate tensile strength 385 1990 
(MPa) 
Elongation at break (%) 10 - 
Density; p (Kg/m?) 2700 7800 
Conductivity; K (W/mC) 90 25 
Expansion coefficient; a  21,1.1 (a 12,6.107© 
(Cc!) 
Specific heat; Cp (J/KgC) 900 500 


interface between the tool and the workpiece. Depending 
on the contact pressure, two values of contact conductance 
were defined: 5000 W/m?C, for the high-pressure condi- 
tions developed during the plunge and the dwell of the 
tool, and 1000 W/m?C, for the low-pressure conditions 
developed during the retraction.*° Thermal properties, 
such as conductivity, specific heat, and thermal expansion 
were also included in the 3D FE models (Table 2). In add- 
ition to the parameters above mentioned, a fraction of 0.9 
was used to define the amount of heat generation asso- 
ciated with plastic deformation. For the tool-workpiece 
contact, it was assumed that 80% of the frictional heat 
was dissipated into the workpiece. For the mechanical 
boundary conditions, a plunge rate of 0.0417 m/s, and a 
rotational speed of 600 rpm were prescribed for the tool. 
The plunge depth was set equal to 0.5 mm. To represent 
the fixture and clamping system, the bottom and lateral 
surfaces of the workpiece were constrained; so that no 
lateral movement was allowed. 


Experimental set up 


Experimental spot welding was performed by plunging a 
pinless rotating tool into the workpiece, until reaching a 
predefined depth (position control), and then maintaining 


the tool at the same position, for a dwell period of 60 s. 
The tool used in these experiments was made of H13 
steel. The material properties of the tool and the work- 
piece are listed in Table 2. In order to compare numerical 
results to those obtained in experiments, the same rota- 
tional speed (600 rpm) and plunge depth (0.5 mm) of 
the numerical simulations were used in the experimental 
tests. During the welding operation, the temperature evo- 
lution was registered using a thermographic camera, fol- 
lowing the procedures explained by Andrade et al.?’ 
After welding, cross-sections were prepared by cutting 
across the centre of the welds to identify the different 
Process Affected Zones (PAZ). Standard techniques 
were used in preparing the metallographic samples for 
optical microscopy. 


Numerical results and discussion 


Temperature distribution 


Figure 3 compares the temperature distribution obtained 
with the different friction models, during the plunge, the 
dwell, and the retraction stages. By examining the differ- 
ent graphs, it can be observed that the size of the high- 
temperature zone and its evolution along with the 
welding time are changed by varying the friction model. 
When using the CLB model, the temperature increases 
during the plunging and spreads throughout the contact 
interface to reach a maximum value of 300°C at the 
start of the dwell period. However, before the end of the 
dwell period, a significant decrease in the maximum tem- 
perature occurs. On the other hand, when using the 
MDF-CLB model or the COF(T) based model, the 
highest temperature region does not spread to the entire 
contact interface and the peak temperature continuously 
increases until the end of the dwelling stage. Moreover, 
compared to the CLB model, it can be clearly seen that 
lower peak temperatures were computed, i.e. 250 °C 
using the MDF-CLB model and 225°C using the 
COF(T) based model. At the beginning of the retraction, 
independently of the friction model, temperatures 
exhibit a gradual decrease, as a result of the convective 
heat transfer between the plate and the surrounding area. 

Figure 4 displays a close-up view of the temperature 
distribution calculated at the end of the dwelling stage 
for the three friction models. The experimental weld 
micrograph was given for comparison purposes. 
Analysing the figure, it is obvious to observe that the 
highest temperature region was concentrated at the 
contact interface between the tool and the workpiece. 
For better lucidity, the contour of this region was overlaid 
with the PAZ (dashed line). Its shape and size coincide 
well with the PAZ region shown in the micrograph, parti- 
cularly for the MDF-CLB and COF(T) based models. 
Both predicted temperature distribution and experimental 
cross-section weld present a W-shape contour. Similar 
findings concerning the temperature distribution and 
weld micrograph were reported by Yang et al.?° The tem- 
perature distribution in the weld zone is a result of the 
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CLB model 


MDF-CLB model 


COF(T) based model 


During the plunge 


Beginning of the dwell 


End of the dwell 


During the retraction 


Figure 3. Comparison of simulated temperature distributions in the weld cross-section using different friction models during the 


different stages of FSSW process. 
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150.00 
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0.00 
MDF-CLB model 


Experimental micrograph 


2mm 


CLB model 


Figure 4. Comparison between simulated temperature distribution by the end of the dwelling stage and experimental micrograph in 


the weld cross-section. 


combined effect of heat generation and dissipation during 
FSSW. The heat is generated through two different 
thermo-mechanical mechanisms, namely the friction 
work and the plastic deformation work. According to 
the numerical results, the highest temperature zone is 
located right below the tool, but away from the weld 
centre. This result could be attributed to the use of a 
pinless tool with a flat shoulder. The absence of the pin 


stirring action results in reduced material flow in the 
weld centre, which leads to a decrease in the heat gener- 
ated from plastic deformation work. 

To illustrate this interpretation, Figure 5 displays the 
equivalent plastic strain distribution in the radial direction 
at the top surface of the workpiece for all the numerical 
simulations. A simple visual inspection of the curves 
reveals that the shape of the equivalent plastic strain 
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distribution is similar when using different friction 
models. The highest deformation values are always regis- 
tered near the tool edge and become extremely low toward 
the centre. However, a significant difference in the peak 
equivalent plastic strain is observed for the CLB model, 
which may explain the difference in the maximum tem- 
perature reached with this model. 

Previous results showed that the temperature evolution 
during the simulation of FSSW was deeply influenced by 
the friction model. However, the mechanisms behind this 
evolution require further investigation. Therefore, in sub- 
sequent sections, we will focus on the study of the contri- 
bution of the energies dissipated by friction and plastic 
deformation on the heat generation. 


Maximum temperature evolution with dissipated 
energy 


As mentioned earlier, in FSSW the temperature evolution 
is highly sensitive to the amount of heat generated through 
two distinct mechanisms: the friction work between the 
contact surfaces (sliding condition) and the plastic 
deformation of the workpiece material (sticking condi- 
tion).’* In this section, the contact developed during 
welding will be assessed by evaluating the friction and 
plastic dissipated energies. Figure 6 shows the numerical 
results of the peak temperature evolution, with the dissi- 
pated energies, for the different friction models. The dif- 
ferent stages of the FSSW process are identified in each 
graph. It can be seen that the predicted maximum tempera- 
ture and energy dissipation are dependent on the used fric- 
tion model. Thus, the effect of each friction model on the 
numerical results is discussed in the following. 


The conventional Coulomb model (CLB). When the CLB 
model was used (Figure 6(a)), during the plunging 
stage, the temperature reached a maximum value of 
200 °C, when the plastic and frictional dissipated energies 
attained the maximum values of 75 J and 5 J, respectively. 
At the beginning of the dwell, the peak temperature 
increased continuously with the plastic dissipated 
energy, meanwhile, the frictional dissipated energy 
remained in very low values. This feature reflects that 
the full sticking condition was developed at the contact 
interface. After a lapse of time, the temperature remained 
almost constant, despite the increase of the plastic dissi- 
pated energy. This result is attributed to the heat conduc- 
tion transfer throughout the workpiece. In fact, with the 
welding time, the generated heat was spread toward the 
bottom surface of the workpiece. At the end of the dwell- 
ing stage, the plastic dissipated energy tended to stabilize, 
while a slight increase in the frictional dissipated energy 
took place. This indicated the beginning of sliding at the 
contact interface (Teri, = 0.301). However, the value of 
the peak temperature started to decrease. This result 
allows deducing that the temperature evolution during 
the process was highly dependent on the plastic dissipated 
energy. 


The modified Coulomb model (MDF-CLB). Addressing the 
results relative to the MDF-CLB model (Figure 6(b)), 
from a general point of view, it is possible to conclude 
that an important amount of frictional energy is dissipated 
when compared to the CLB model. This result indicates 
the presence of the sliding condition in some regions at 
the contact interface due to the introduction of the limit 
shear stress value of the base material. During the plun- 
ging period, the amount of the plastic dissipated energy 
was still higher than that of the frictional dissipated 
energy. It indicated the predominance of the sticking con- 
dition with the downward motion of the tool. At this stage, 
110°C maximum temperature was reached, when the 
attained plastic and frictional dissipated energies were 
58 J and 12 J, respectively. 

It can be noted that, compared to the CLB model, the 
small decrease in the plastic dissipated energy induced a 
sharp decrease in the maximum temperature, regardless 
of the increase of the frictional dissipated energy. When 
the plastic heat generation was localized near the top 
surface of the workpiece, the plastic deformation is 
likely the determinant factor for the temperature evolu- 
tion, mainly during the plunging. Once the dwelling 
stage starts, the temperature increases slowly. It is clear 
that the plastic dissipated energy was steady, meanwhile, 
the frictional dissipated energy was increasing. The 
maximum temperature value reached at the end of the 
plunge stage is still insufficient for softening the work- 
piece. Therefore, with the lack of tight contact between 
the tool and the workpiece, the plastic deformation will 
be trivial and the sliding condition will prevail over the 
sticking one. After some time, when the temperature 
exceeded the value of 150°C (softening temperature), 
the material of the workpiece has experienced excessive 
plastic deformation (Figure 2), which induced a signifi- 
cant rise in the plastic dissipated energy. Thereby, both 
energy components have started increasing at the same 
rate, which indicates that mixed sliding-sticking condi- 
tions were developed. 


The temperature-dependent friction coefficient-based model 
(COF (T)). Figure 6(c) shows the evolution of temperature 
and frictional and plastic dissipated energies for the 
COF(T) based model. It is noted the important contribu- 
tions of both dissipated energies in the heat generation. 
It also indicates the presence of mixed sliding-sticking 
contact conditions. However, compared to the 
MDF-CLB model, a higher temperature was reached at 
the end of the plunge (165 °C), when the plastic and fric- 
tional dissipated energies were 69 J and 9 J, respectively. 
It is also obvious to note that the increase in the tempera- 
ture was caused by the rise in the energy dissipated from 
the plastic deformation work. During the dwelling period, 
a sharp decrease in the rate of plastic dissipated energy 
was observed, contrary to the frictional dissipated 
energy which was increased most rapidly. Similar to the 
two previous models, the reduction in the tool pressure 
after the transition from the plunging to the dwelling 
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Figure 5. Plots of equivalent plastic strain versus distance away from the tool centre for the different friction models. 
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Figure 6. Temperature evolution with dissipated energy for different friction models: (a) CLB model; (b) MDF-CLB model; (c) 
COF(T) model. 
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stage induced a decrease in the plastic deformation along 
the thickness direction. Contrary to the MDF-CLB model, 
the plastic dissipated energy in this case, continued to 
increase from the beginning of the dwelling stage. 
According to Figure 2, the temperature value of 165 °C, 
reached at the end of the plunge stage was high enough 
to cause softening within the weld zone. 

The maximum stabilized temperature, predicted using 
the COF (T) based model, was lowest than that forecasted 
using the other models. In fact, when using the 
temperature-dependent-COF based model, different 
values of critical shear stress were considered according 
to the temperature reached at the contact interface. This 
fact allows the sliding condition to overcome the sticking 
one, which is reflected by the high amount of the frictional 
dissipated energy during the process. Moreover, as pre- 
sented in Table 2, the values of COF given by 
Meyghani et al.?* were lower than 0.3 (introduced in the 
CLB and the MDF-CLB models), which leads to an 
earlier transition from the sticking to the sliding state. 
These findings for the COF (T) based model are consistent 
with Su et al.** works. It was stated, based on an analytical 
study for the FSSW process, that almost the whole 
mechanical energy was converted into viscous heat. 
Schmidt & Hattel?? drew similar conclusions, assuming 
that the sticking condition predominated during welding. 
Thus, the heat generation at the tool-workpiece interface 
was primarily caused by plastic dissipation energy. 
Similar conclusions, regarding the heat generation propor- 
tions, were reported by Li et al.>! The authors claimed that 
the heat generation during FSW depends on the intense 
plastic deformation of the thin shear layer at the tool- 
workpiece interface. 

It should also be noted that, for all the friction models, 
a steady-state of the dissipated energies during the retrac- 
tion stage is explained by the fact that the energy extracted 
from ABAQUS® is a cumulative result (history output). 
Therefore, when there is no energy source, the dissipated 
energies remain constant. 


Numerical versus experimental results 


The analysis performed in the previous item showed the 
contact conditions prescribed according to the different 
friction models, as well as the temperature evolution 
during welding. The predicted temperature evolution 
was compared to the experimental results, to examine 
the accuracy of each friction model. Figure 7 shows the 
predicted and experimental maximum temperature evolu- 
tion. One can observe that the numerical simulation per- 
formed using the COF (T) based model is in good 
agreement with the experimental measurements. 
Compared with experimental results, the thermal cycle 
obtained using the CLB model displays a maximal 
average relative error of 50% (about 100 °C). This result 
disproves the assumption of 0.3 as a COF value for the 
present process conditions. According to Dixit et al.°? 
any friction model can be linearized in the form of a 
CLB model for a very narrow load range; but further 
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Figure 7. Comparison of the maximum temperature evolution 
using different friction models with experimental results. 


than a certain load, the COF decreases along with the 
normal pressure. This statement may explain the very 
large differences in the predicted temperature between 
the CLB model and the experimental results. Therefore, 
one should be aware in the selection of the COF when 
using the CLB model. 

For the MDF-CLB model, the difference between the 
numerical and the experimental results decreases, reach- 
ing an average relative error of about 25% during the 
dwell period. It is important to stress that, during the 
plunge period until the beginning of the dwelling stage, 
the model underestimates the welding temperature. 
However, during the dwell stage, the model overestimates 
the temperature, once the plastic dissipated energy starts 
to increase. These results indicate the weakness of 
MDF-CLB models, mainly when the material exhibits 
large plastic deformation. It is also important to mention 
that the yield shear stress (t,je¢) of the material during 
welding depends on the temperature. However, the evolu- 
tion of t,;eq with the temperature is not considered when 
defining the MDF-CLB model under ABAQUSC©. This 
was one of the software limitations. Despite the limitation 
of the conventional and the modified versions of the 
Coulomb model in the rational description of friction in 
the numerical simulation of FSSW process, they are the 
most used in literature. Contrary, we showed that the 
COF (T) based model predicts high accurate results com- 
pared to those obtained using CLB or MDF-CLB models 
in the simulation of FSSW process. From this comparison, 
it can be deduced that the friction coefficient evolution 
with temperature (Table 2) reflects the real effect of the 
material softening on the contact state during FSSW. 
This result was endorsed by Moore & Tegart*® who 
claimed that the coefficient of friction decreases with the 
yield pressure. 

In order to improve the accuracy of the CLB and 
MDF-CLB friction models in the prediction of the 
welding temperatures, a simple approach was experi- 
enced. More precisely, the constant COF, for both 
models, was adjusted to a lower value (u=0.23). The 
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Figure 8. Numerical results obtained with CLB model and MDF-CLB model when using a COF = 0.23: (a) The maximum temperature 


evolution; (b) frictional and plastic dissipation energies. 


new results of maximum temperature evolution are 
plotted in Figure 8(a). The comparison between the com- 
puted and experimentally obtained results shows a clear 
improvement in the predictive capability of both friction 
models. 

Figure 8 displays the results of the frictional and plastic 
dissipated energies using the friction coefficient of 0.23. 
Compared to the results presented in Figure 6(a) and 
(b), a significant decrease in the amount of plastic dissi- 
pated energy is detected. Meanwhile, the frictional dissi- 
pated energy exhibits a significant increase. This result 
reflects the earlier occurrence of the sliding condition 
when reducing the COF value. According to Mekote 
et al.*! the apparent friction coefficient measured using 
tribometers is often overestimated, because they did not 
consider the intense contact pressure and the plastic 
deformation that were taking place during welding. 
Moreover, depending on the relative motion of the 
contact surfaces, two types of friction may occur: 
kinetic friction associated with a linear motion, and 
rolling friction associated with rotational motion. The 
coefficient of rolling friction is generally smaller than 
the kinetic friction coefficient.** Therefore, to improve 
accuracy in the FSSW simulation, the contact model 
should involve the component of the rolling friction; 
thus a reduced value of the COF is to be considered. 
Nevertheless, the effect of temperature on the friction 
coefficient should also be taken into account. 
Consequently, an accurate prediction and a better com- 
parison with experimental results could be achieved. 


Conclusions 


This work puts forward a thorough analysis of the effect 
of three friction models on the heat generation sources 
(i.e. energy dissipation by friction and plastic deform- 
ation) and the temperature evolution in the FSSW 
process simulation. A 3D thermo-mechanical FE analysis, 
using three different friction models, namely, the conven- 
tional Coulomb model (CLB), the modified Coulomb 


model (MDF-CLB), and the temperature-dependent fric- 
tion coefficient-based model (COF(T)), was performed 
for the AA6082-T6 aluminium alloy. For each friction 
model, the amount of energies dissipated by friction and 
by the plastic deformation along with their contribution 
to the maximum temperature evolution were calculated. 
Then, the comparison between prediction of the models 
and the experiments was established. Based on this ana- 
lysis, the following conclusions can be drawn: 


e For the AA6082-T6 aluminium alloy, which exhibits 
large plastic deformation during welding, the plastic 
dissipated energy was the main source of heat gener- 
ation and subsequently the determining factor in the 
evolution of temperature during the process. 

e During the FSSW process, the predominance of a 
contact condition, whether sliding or sticking, affects 
the amount of heat generation from both friction 
work and plastic deformation. Therefore, the friction 
models parameters, mainly, the friction coefficient 
(COF) and the material yield shear stress (tyie), 
which determine the transition from the sticking to 
sliding condition, need to be fine-tuned. 

e When using the CLB and the MDF-CLB friction 
models, the constant values of COF and Tyjeiq should 
consider the high contact pressure and the severe mate- 
rial plastic deformation. The experimentally measured 
COFs are often overestimated, and reduced values 
need to be introduced for more reliable numerical 
simulations. 

e The COF(T) based model, which considers the material 
softening at high temperature, was found to be the most 
suitable friction model to predict the real temperature 
evolution during the FSSW of AA6082-T6 aluminium 
alloy. 


These findings clearly demonstrated that an efficient 
numerical analysis of FSSW may be achieved when con- 
sidering the material behaviour in the selection of the fric- 
tion models parameters. This approach could eventually 
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lead to an efficient contact modelling strategy for the 
numerical simulation of the FSSW and the FSW 
processes. 
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